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Abstract 

Using appropriate harmonics, we study the future asymptotic behav- 
ior of massless scalar fields on a class of cosmological vacuum spacetimes. 
The spatial manifold is assumed to be a circle bundle over a higher genus 
surface with a locally homogeneous metric. Such a manifold corresponds 
to the SL(2, R)-geometry (Bianchi VIII type) or the H 2 x R-geometry 
(Bianchi III type). After a technical preparation including an introduction 
of suitable harmonics for the circle-fibered Bianchi VIII to separate vari- 
ables, we derive systems of ordinary differential equations for the scalar 
field. We present future asymptotic solutions for these equations in a 
special case, and find that there is a close similarity with those on the 
circle-fibered Bianchi III spacetime. We discuss implications of this simi- 
larity, especially to (gravitational) linear perturbations. We also point out 
that this similarity can be explained by the fiber term dominated behavior 
of the two models. 

1 Introduction 

Linear perturbation analysis of Einstein's equation has fundamental importance 
in general relativity. It expands the significance of an exact solution that usu- 
ally has a large symmetry by providing additional properties of the solution, 
especially those about stability or instability when the solution deviates from 
the symmetric configuration. In particular, detailed properties of the perturba- 
tions are expected to be served as useful pieces of information with which global 
nonlinear analyses, like those for the global existence problem [21 EH an d the 
cosmic censorship conjecture [181 114) or conformal dynamics by Einstein's flow 
00], can develop. 

Linear perturbation analysis can be however very difficult and often requires 
a formidable effort to carry out for many of the solutions we are interested 
in. Among them are the spatially homogeneous solutions [HI EH which are not 
isotropic in any limit keeping the homogeneity. Linear stability properties are 
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not known for most of those solutions. Recently, however, a complete pertur- 
bation analysis for a vacuum Bianchi III solution, which belongs to this class 
of anisotropic solutions, has been carried out [213 E2] after an enormous effort. 
In connection with this, we in this paper propose an indirect approach that 
extends the Bianchi III result to another model. More precisely, we make use 
of an analogy to derive perturbation properties for a "Bianchi VIII" solution. 

The setting is as follows. The spacetime manifold we consider is the direct 
product of a spatial manifold M and time R, for which M is a circle bundle 
over a higher genus surface £ s with genus g > 1 . The background spatial metric 
is assumed to be locally homogeneous. In this situation, there are two kinds 
of possible locally homogeneous metrics on M depending upon whether M is 
a nontrivial bundle or a trivial bundle. If the bundle is nontrivial the metric 
becomes the Bianchi VIII type (Thurston's SL(2, R) [H|[2Hj), while if it is trivial 
it becomes the Bianchi III type (Thurston's H 2 x R). 

The key fact is that the two kinds of vacuum spacetimes behave the "same" 
way in the future asymptotics in a special case; the evolutions of the fibers 
become asymptotically the same, and the evolutions of the base also become 
asymptotically the same. We will call this coincidence of the two backgrounds 
the background asymptotic degeneracy. (This "degeneracy" makes sense only 
when considering a class of fibered spatial manifolds with the same fiber and 
base manifold.) With this kind of coincidence, it is of great interest to compare 
test fields like scalar fields on these spacetimes and see if there is a substantial 
difference in the asymptotic behaviors of the fields. That is, if we find that there 
is no substantial difference in those behaviors — we will call this coincidence the 
scalar field asymptotic degeneracy, then this may be regarded as an evidence that 
other kinds of fields like the linear perturbations on the two kinds of spacetimes 
will also show the same asymptotic behavior or at least a very similar behavior 
to each other. 

Based on this idea, in this paper we study the massless scalar field equation 
on the Bianchi VIII spacetime in detail. In particular, we separate the field 
equation using appropriate harmonics and find asymptotic solutions. We then 
make a comparison with the corresponding Bianchi III system, and will confirm 
the scalar field asymptotic degeneracy between the two systems. As mentioned, 
since the basic properties of the perturbations of the Bianchi III solution are 
already known, the comparison provides us with a conjecture about the Bianchi 
VIII perturbations. 

One of the unfamiliar techniques needed to carry out our analysis may be 
the harmonics for the circlc-fibered Bianchi VIII manifold. Since there does not 
seem to exist a prior work on this subject (see however |4*llll| for related work), 
we start with introducing appropriate ones for this manifold in §§ 12131 Then, we 
separate the field equation using those harmonics and find asymptotic solutions 
of the reduced equations in §5 14151 We then compare the results with those of 
the Bianchi III system in §03 The final section is devoted to conclusions. 

2 SL(2,R) actions 

In this section we introduce some basic objects associated with the Bianchi VIII 
manifolds as a preliminaxy for developing the harmonics. 

Let G — SL(2,R) be the universal covering group of SL(2,R), the real 
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special linear group of rank 2. This group is also known as the Bianchi VIII 
group. Our spatial manifold M is the quotient T\G of G by a discrete subgroup 
r C G acting freely from the left. M is assumed to be closed (— compact 
without boundary). In such a case, the manifold M is known to become a (non- 
trivial) circle bundle over a higher genus hyperbolic surface Sg. 1 (g > 2 is the 
genus.) 

For given element 

9=( a c b d ) eSL(2,R), (1) 

it acts on the space of direct product fi + x T of the upper half plane 

n + = {( € C|ImC > 0} (2) 

and the circle T = R/2ttZ, 

j:!1 + xT^!!|XT, (3) 



as 

s-(C,*) = (^(0,*-**(0), (4) 

where (£, z) — (x + iy, z) E f2 + x T, and 

^(0 = ^±^, ^(C) = 2arg(cC + rf). (5) 
In fact, it is a direct computation to confirm the homomorphism 

(gg r ) • (C,*) = 9 ■ (</ • (£*)), (6) 

where 173' is understood as the usual matrix product. Since this natural action 
is simply transitive, i.e., for arbitrary x, x' G Q + x T there exists unique g € 
SL(2, R) such that g ■ x = x', we can identify the space f2+ x T with the group 
SL(2, R) by associating g S SL(2,R) with its action on the origin Xo = (0+li,0); 
<7 ~ g • Xo. If we take the universal cover of 17+ x T, it apparently becomes 
f2+ xR, which can be also identified with the space R + = {(x, y, z) £ R 3 \y > 0}, 
therefore 

G = SL(2,R)~rj+ xR~R + . (7) 

This identification is understood throughout this paper. The following explicit 
formula for the correspondence G ~ R+ is sometimes useful; 



Vv V 



(x,y,z). (8) 



To find generators of G acting on R, , it is convenient to introduce the 
following one-parameter subgroups of SL(2, R); 

1 s \ ( e*/ 2 \ / cos 6>/2 sin 6/2 \ , n . 

1 ) ' at - { e^*/ 2 J ' " e ~ I -sin^/2 cos0/2 J ' (9) 



^^More precisely, M can be a Seifeld fiber space over a hyperbolic orbifold 1191 . but in this 
paper we assume that M is one of usual fiber bundles as a typical case. 
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where £ [0, 4tt), t € R, and seR. Then, any element 5 6 SL(2, R) is uniquely 
expressed in the form (|U, Proposition 1.3, Chap V) 

g = u e a t n s . (10) 

This decomposition is known as the Iwasawa decomposition. Let £1, £2, and £3 
be the generators of n s , a t , and ug, respectively; so, e.g., for x e R 3 ^ , 

6/(x) = ^/K.x)| s=0 = A /(x) . (11) 

Similarly, we can compute the others, and find 

f 9 
6 ~ "a - ' 

Ox 

d d 
1 1 2 2 ,\ ^ 9 d 

These generators satisfy the following algebra; 

= a, = -a + 6, [6,6] = -&. (13) 

Although this is not of the canonical form that is usually assumed as the algebra 
of Bianchi VIII, it is easy to see that the linear combinations 



a = ^(-^1+6), £2=6, £3 = ^(56+6), (14) 
satisfy 

[&,&]= C K „& (15) 

with the structure constants C K jj = —C k ji for which the nonzero independent 
components are given by 

C\3 = C 2 31 = 1, C\ 2 = -1. (16) 

This coincides with the canonical choice of structure constants of Bianchi VIII 
(e.g., 0). 

The generators 6 above are regarded as generating the left action of G on 
G ~ R 1 ^. Another important kind of objects are generators of the right action, 
which we denote as \i ■ They are defined as differential operators such that they 
commute with all the left generators 6 

= 0, (I,J = 1-3), (17) 

and possess the same structure constants ()16[). They are given by 

( . d d \ . d 

Xi=V\ smz ^ cosz "5~ ~ smz "5-' 

y ox oy J oz 



d 

oz 



(18) 
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The commutativity (|17|l means that xi are invariant under the left action of 
G as seen from the vanishing of the Lie derivatives, C^XJ — 0. We therefore 
call xi the invariant vectors. The dual one-forms of these vectors, the invariant 
one-forms, can be expressed as 

a 1 = — (sin zdx — cos zdy), 
V 

a 2 = — (cos zdx + sin zdy) , (19) 
V 

o dx 

er = 1- dz. 

V 

These one- forms are invariant in the sense C^ I a J = 0, and satisfy the following 
Maurer-Cartan relation 



da 1 = -\c I JK a J Au K , 



(20) 



with respect to the same structure constants (|16|) . Now, a homogeneous metric 
on G can be expressed as q = qi.jcr 1 cr J with qu being constants, since we 
apparently have C^, q — for this metric. The generators of the left actions 
are Killing vectors for these homogeneous metrics. When the metric is of the 
form q = ^((cr 1 ) 2 + (cr 2 ) 2 ) + 93 (cr 3 ) 2 , where qi and 93 are constants, this metric 
possesses a fourth Killing vector 

& = 3T- ( 21 ) 

and the metric is said to be locally rotationally symmetric (LRS). Note that this 
additional Killing vector for an LRS metric coincides with one of the invariant 
vectors, X3 5 

£4 = X3- (22) 

An importance of the invariant vectors Xi is that they generate the regular 
representation (T, L 2 (T\G)) of G on the Hilbert space L 2 (r\G) of all square- 
integrable functions on T\G. By regular representation, we mean the "right 
regular representation." 2 It is apparent that xi generate the regular represen- 
tation on L 2 (G), since xi are generators of the right action (of G) on G. Note 
that the invariant vectors xi are naturally well defined on the quotient T\G, as 
well, since xi are by definition invariant under the action of V C G. (Therefore 
for the covering map 7r : G — > L\G, the induced vectors tt*xi on T\G are well 
defined, which however we simply denote as xi-) Moreover, since the right and 
left actions commute each other the right action of G on the left quotient T\G 
is also well defined. Therefore the regular representation (T, L 2 (T\G)) is well 
defined and xi are its generators. 

With the aid of Eq.(|SJ), it is easy to explicitly compute the right action x • g 
of g 6 SL(2,R) on x e R^. (Compute the usual matrix product of the left 
hand side of Eq.JHJ) and g, and then read off the components in R + , using 

2 For given g £ G, we can define the map T g : /(x) — + /(x ■ g), where x £ G ~ ~R? + . The 
right action x • g is naturally defined viewing x as the corresponding element in G. It is easy 
to confirm that the map T : g — » T g is a homomorphism; T g T g i = T„ g i. This homomorphism 
(T, L 2 (G)) is called the (right) regular representation of G on L 2 (G).^l) 
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the correspondence ©.) The most important action is that of the compact 
subgroup U0, which is found to cause the translations (x,y,z) — > (x,y,z + 9). 
Its generator is therefore found to coincide with X3> which we call the fiber 
generator. 

Ea. (|22fl is saying that if the metric is LRS, the fiber generator X3 is also a 
Killing vector. Therefore an LRS metric of Bianchi VIII type is U(l)- symmetric 
(as well as being locally homogeneous), since there exists isometries along the 
circle (~ U(l)) fibers. Beware that non-LRS Bianchi VIII metrics are not U(l)- 
symmetric in any sense despite that they are circle-fibered and locally homoge- 
neous. This is apparent from the fact that none of the generic Killing vectors 

(/ = 1 ~ 3) are well defined on the quotient manifold T\G, since £/ are not 
commutative with each other. (Remember that the action of T is generated by 
£j. This noncommutativity therefore implies the noncommutativity between T 
and £/, which in turn implies £j are not well defined on T\G.) 

The regular representation (T, L 2 (T\G)) is, as shown in the next section, 
unitary, but not irreducible. As well known (e.g., [2]), harmonics are basis 
vectors of irreducible components of a regular representation. We are hence 
interested in irreducible components in (T, L 2 (T\G)), which are the subject of 
the next section. 



3 The harmonics 

The most important entity to consider an irreducible unitary representation 
is the Casimir operator, which is an operator which commutes with all the 
generators of the representation, since from Schur's lemma, such an operator 
must be a constant when acting on an irreducible space. For the group G, it is 
given by 

□ ^(Xi) 2 + (X2) 2 -(X 3 ) 2 . (23) 

In fact, it is a direct computation to check the condition [□, X/] — for all 
I = 1 ~ 3. In particular, it commutes with the fiber generator \3', 

P,Xs]=0. (24) 

It is therefore possible to simultaneously diagonalize these two operators; 

□<£m,A = ~A0 m A) 
X3<Pm,A = im0 m ,A- 

We call A the Casimir eigenvalue, m the fiber eigenvalue. 

We assume that the closed manifold T\G coincides with a compactification 
of SL(2, R) (rather than G = SL(2, R)). This is equivalent to assuming that the 
discrete subgroup T has the subgroup 2nZ generated by the action (x, y, z) — > 
{x,y,z + 2ir). Any two points (x,y,z) and {x,y,z + 2tt) should therefore be 
identified, which forces 

m e Z. (26) 

If we considered a p-fold covering of the manifold, we would instead have 
m G Z/p = {n/p\n 6 Z}. For simplicity, however, we do not consider this 
generalization in this paper. 
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We assume that each <p m ,A is normalized to unity; 

||0m,A|| 2 = (0m, A, 4>m,h) = 1, (27) 

with respect to the inner product (/i, fy) in L 2 (r\G) defined by 



(/i,/ 2 )= / fifidfio. (28) 
Jr\G 

Here, /* is the complex conjugate of /. The measure defined by 

dfJ, = a 1 A a 2 A a 3 (29) 

is the natural left-invariant measure of G, therefore the above integral on the 
quotient T\G is well defined. Since the measure is also right-invariant (i.e., it is 
"unimodular" ) , the above inner product is invariant under the right action of 
G, i.e., for g e G, 

/ h(x-g)fZ{x>g)dno= I /i(x)/ 2 *(x)d Mo . (30) 
Jr\a Jr\c 

This implies that the regular representation (T, L 2 (T\G)) is certainly unitary. 

Our main purpose here is to find the harmonics belonging to given A, i.e., to 
find all the basis functions {<f) m > \} m > of an irreducible space labeled by A. (The 
number of the irreducible spaces belonging to the same A, the "multiplicity," 
can be more than one. So, A is not a complete label to specify the irreducible 
space. More about the multiplicity will be mentioned at the end of this section.) 
This task of finding the harmonics can be done starting from one function <j) mt A 
satisfying Eqs.J2HJ), as shown below. 

First, let us define 

We have the following properties about adjointness: 
Lemma 1 In the Hilbert space L 2 (F\G), it holds 

4 = -A 2 , A\ = A 3 . (32) 

Proof. We first show that the invariant operators xi (I = 1 ~ 3) are 
anti-selfadjoint; \\ = — Xi- Actually, if it is the case, A\ = ^g (xi — iX2^ = 

^(-Xi - iX2) = -A2, and A\ = -{ixiV = -«X3 = A 3 - To show the anti- 
self adjointness of xii note the identity 

{Xih,h)=Ti-{h,Xih), (33) 
where we want to show the integral 



1/ = / Xi(fifZ)dlH> (34) 
Jr\G 

vanishes for all I = 1 ~ 3. Actually, we have 

1i= [ Xitfi&V A a 2 A a 3 = \ [ d{h!^UKO J Aa K ), (35) 
Jr\G z Jr\G 
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which is confirmed using the identity 

df = {xifV + (X2fV + (X3/)a 3 , (36) 

and the relation 

d{e IJK a J A a K ) = 0, (37) 

where sijk is the unit skew symmetric symbol; cjjk — C[ijk], £123 = 1- Then, 
from Stokes theorem the quantity 2j does vanish, since the manifold T\G is 
assumed to be closed. ■ 

Next, note that the commutation relations among Ai become 

[A 3 , Ai] =Ai, [A3,A 2 ] = -A 2 , [Ai,A 2 ] = A3. (38) 

This in particular shows that operators A\ and Ai are, respectively, the raising 
and lowering operator. In fact, since 

A 3 AKpm = ([A 3 , A] + A-4 3 )0m = {Ai + AiA 3 )<j> m = (1 + m)Ai<p m , (39) 

A\4> m is an eigenfunction for m' = m + 1; A\4> m oc 4> m +i- Similarly, A24> m oc 
4> m -x- We can therefore assume that 

Ai4> m — a m (f> m+ i, 

Azfym = b m (j) m -i, (40) 
A3(j) m = mifim, 

for appropriate coefficients a m and b m . 

On the other hand, one can easily check the identity 

a 1 a 2 ^1(d-(A3) 2 + a 3 ). (4i) 

Together with the assumption (|4T1)l we find 

a m -xb m = Z ^ Um-i) 2 +A-ij . (42) 



From Ea. H32Jl . we have 



(AlA24>m,h,4>m,k) = {A2(j} m ,A, ~A2(t) m ,A) 

= -||^m,A|| 2 (43) 

= -\b„ 



'm\ , 



(44) 



where we have used the normalization (|27|) . On the other hand, 

{AlA24>m,A, 4>m,A) = tt m -l& m (^m,A. 0m, a) 

Therefore, using Ea. ll4*2l) . 

\b m \ 2 = -a m - X b m = i Urn - i) 2 + A - i J . (45) 
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In particular, since \b m \ > 0, for any possible m and A, it must hold that 

(m-i) 2 +A-l>0. (46) 

Since an m = mode is always contained in any irreducible components, 3 from 
the above condition A should be nonnegative 

A > 0. (47) 

Also, 

b* m = -Otn-i. (48) 

We can determine the coefficients a m and b m from Eas. i|45[) and (|48|l . For 
convenience, let us define 

s ee ±J\A-~\. (49) 
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(i) Case A > \ 

In this case, it is convenient to choose 



1 n !| ■ x 



(50) 



(ii) Case < A < \ 

In this case, we can choose 



2' 



1/2 
1/2 



(51) 



When A = (s 2 = 1/4), from Eqs.JSTJ we have ao = 60 = 0, reflecting the 
fact that this representation is trivial (see below) and the representation space 
is spanned by only one constant function 4>o.o- The other representations are all 
infinite dimensional. 

The differential representations H4()|l with the coefficients H5()|l or (|51|) above 
provide the heart of the relations used to separate the field equations. 

In the classification of irreducible unitary representations of SL(2,R), there 
are five kinds of series (e.g., |211 \'27\ ). The case (i) above belongs to the class 
called the first principal series, while the case (ii) the complementary series. 
The other series may not generically occur from the representation we are con- 
sidering, for which the fact (|2(j|l holds. 

3 This can be seen from Eq. I26I and the fact that Ai and A.2 raises or lower the fiber 
eigenvalue by 1. That is, Spec.43, the spectra of ,43 for given A, always coincides with the 
whole range of possible m; Spec.43 = Z. (A = is an exceptional case, but m = exists in 
this case, too.) Remark however that this would not be the case if m (i Z/p (p > 1), which 
was as mentioned possible if considering a covering. 
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Note that the recursion relations l|4(J|l imply that <fi m ,A can be constructed 
by successively applying appropriate operators on 4>o,a- Explicitly, the relation 
is given by 



Dm, A 



(V-l^)(«m-2^) • • • K Ui)^o, A (m > 0) 

fei^XC*^) • • • (^-Wo.A (m < 0). 



The function 4>o.A is a z-independent function, since from definition X3</>o.a = 
d<f>Q^\/dz = 0. Moreover, let us observe the fact that when acting on a z- 
independent function f(x,y), the Casimir operator □ degenerates to a two- 
dimensional Laplacian A . (We use hat " for quantities on a surface or z- 
independent functions.) In fact, since we can compute, from the expression l|18(l 
of Xi, 

M > + fe ) W( |l + |\_ 2s ^ + |l (53) 



we have 



where 



□/ = ((Xi) 2 + (X2) 2 - (Xa) 2 )/ 

= ((Xi) 2 + {X2?)f (54) 

= An/, 



The operator A does coincide with the Laplacian, associated with the standard 
two-dimensional hyperbolic metric h = (dx 2 + dy 2 )/y 2 — {cr 1 ) 2 + (cr 2 ) 2 on S g . 

Therefore we can summarize the construction of the harmonics on T\G as 
follows; we first consider a solution of the two-dimensional eigenvalue equation 

A O A = -A0 A (56) 

on the hyperbolic surface (E ff ,/i), and then apply the formula (|52|l assuming 
4>o,A = 0A to obtain an irreducible set of harmonics. Repeat this procedure for 
all independent solutions for given A. Repeat this for all possible values of A, 
then we obtain a complete set of harmonics. 

This construction implies a simple but remarkable fact |27j : the multiplicity 
in an irreducible representation specified by A, occurring in (T, L 2 (r\G)), is 
directly determined by the multiplicity in the corresponding eigenstate of the 
two-dimensional Laplacian A . 

The spectrum of (minus) the Laplacian — A on a hyperbolic closed surface 
has been one of the major subjects in Riemannian geometry 5 . In particular, 
those in the range (0, 1/4) are called small eigenvalues, and in what conditions 
they appear has been one of the central issues. We do not discuss this further 
here, but remark that small eigenvalues appear only when S 9 takes a particular 
"shape" (specified by some "Teichmiiler parameters"), in particular they do not 
necessary exist for an arbitrarily "shape" of £ s . The eigenvalue A = corre- 
sponds to the trivial case, so the multiplicity (in the irreducible representation 
for A = 0) is always one. 
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4 The background solutions 

We assume that the background spacetime metric is (locally) of the form 

g = -N 2 (t)dt 2 + q^a 1 ) 2 + q 2 (t)(a 2 ) 2 + q 3 (t)(a 3 ) 2 . (57) 

No general vacuum solution for this metric is known. However, see |17| for the 
future asymptotic behavior. 

On the other hand, the general vacuum solution for the LRS metric, i.e., the 
case qi(t) = q<i (t), is known. It is a special case of the so-called NUT solutions 
[T5| . The Bianchi VIII NUT solution is given by 

N(t) 2 = U(t)-\ qi (t) = q 2 (t)=t 2 + l 2 , q 3 (t)=Al 2 U(t), (58) 

where 

_ t 2 - I 2 + 2jrf 

u w = — 7^72 — » ( 59 ) 

with I > and \x being constant parameters of solution. (This form of solution is 
obtained as a special case of the metric derived in 0.) Note that the conformal 
metric h = (a 1 ) 2 + (a 2 ) 2 is the hyperbolic metric in the upper half plane; 
h = (dx 2 + dy 2 )/y 2 , thus the conformal factor qi(= q 2 ) is the scale factor for 
the base surface, while q 3 is the one for the circle fibers. To describe the future 
asymptotic behavior of these scale factors, let us introduce a proper time r for 
this LRS solution by 



Ndt = I \J- lj2 dt 
*-A*logt + 0(i). 



(60) 



This implies 



r + /ilogt + 0(i) 

T + /ilog(r + M logt + 0(i)) + 0(i) (61) 

log T 

r + /ilogr + 0(- fc ' 



T 

Then, it is easy to confirm the following expressions. 

\2 



?i(t) =<&(t) = (r + n\ogTf +0(logr), 



T 



In particular, note that the fiber length 

,-2tt 



L= [ f ^dz = 2w^ = 2Tr(2l + 0(-)) (63) 

ifiber JO T 

approaches constant = lmv^oo L = Anl. Therefore the parameter I is 
(47r)~ 1 times the fiber length at the future infinity. 
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5 Future asymptotics of scalar fields 

Let us consider the massless scalar field equation 4 

g ah V a V 6 * = 0, (64) 
where V a is the covariant derivative operator associated with a spacetimc metric 

gab- 

For our shift-free, Bianchi VIII metric 1)57(1 . 

^^={^M^ gN ~ 2 ^) +Aq ) % (65) 



where y/—g = \J — det g a f, — N ^Jqiq^B an d A 9 is the Laplacian with respect to 
the spatial metric q a t, which can be expressed using the invariant operators as 



A ° = H < ?7 1 (X/) 2 - (66) 
i=i 

It is convenient to write the Laplacian A g in terms of Ai and divide it into 

ic 

A ? = A(°) + A« (67) 



two parts, the homogeneous part and the inhomogeneous part Aq\ 



where 

A< 0) = l(q^+q^)(A 1 A 2 +A 2 A 1 ) - q^(A 3 ) 2 , 

2 (68) 
Af^-(q^-q^)((A 1 ) 2 + (A 2 ) 2 ). 

The "homogeneous" part does not change the index m when it acts on a single 
mode function <j) m = 4> m ,M while the "inhomogeneous" part does (and gives rise 
to inhomogeneous terms in field equations). In fact, we find 

A*°Vm = Q^iT 1 + <l2 l )( a ™-ib m + a m b m+1 ) - q^m 2 ^ <j) m 

= -K m (t)ct> m , (69) 

where 

K m (t) = ife- 1 + q^)(m 2 + A) + q^m 2 . (70) 

Here, we have used Eq. l|45|l to deform the first equation of (|69|l . 

We expand the field component W = belonging to an irreducible space 
specified by A as 

$(t,x) = ^^ m (^ m (x). (71) 



*In this section we employ the abstract index notation 1291 and use leading Latin letters 
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Then, the field equation (|64[) reduces to the following equations: 



N d ( (qiq 2 q 3 ) 1/2 



+ N 2 K m (t)^ m =I m {t), (72) 



^.(*) = IT? ( — 7— ) + — (75) 



(<zig 2 <z 3 ) 1/2 dt V ^ dt 

where 

I m (t) = -yCffl -92 )( a m-2am-ll/'rn-2 + 6m+2fem+l'0m+2)- (73) 

These equations form two sets of infinitely many simultaneous ODEs — one for 
{V>m}m=odd and one for {ip m } m=cvcn . The term I m , which contains ^±2, works 
as an inhomogeneous term if we view the above single equation as a dynamical 
equation for -0 m . 

The LRS case. Let us consider the LRS background case (|58|l . in which, 
due to the vanishing of the inhomogeneous term, each mode equation (|72(l be- 
comes independent. The equation becomes of the form 

/ • 

tpm + j4> m + Z rn %l) m = 0, (74) 

where 

■t 2 + l 2 \ 2 . m 2 +A 
and 

f(t) = t 2 -l 2 + 2fJ,t. (76) 

Here, dot (') stands for d/dt. We are interested in the future (t — > +00) asymp- 
totic solution. As emphasized in [22], to this it is necessary to transform the 
time variable to one that is suitable for the analysis. Let us define the new time 
variable T by 

f" = a{t) > °' (7?) 
using a positive function <r(t), which is to be determined. We transform the 
unknown function ip m to another function X m so that the new equation in 
terms of X m , dX m /dT, and d 2 X m /dT 2 has a vanishing dX m /dT term. It is 
easy to do this for given a(t); The transformation is given by 

ijj m = a(t)X m , (78) 

where 

a(t) = a- 1 ' 2 ]- 1 ' 2 - (79) 
The resulting equation is given by 

^L+WX m = 0, (80) 

where 

\ I n> f rv \ 

(81) 

(Again, dot (') stands for d/dt, not d/dT.) If the function W(T) in this equation 
approaches a constant C ^ as T — > 00 (t — > 00), then it may be natural to 
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expect that the equation (|80|) has fundamental solutions approaching e ±l ^ T 

in case C > 0, or e ± ^\^~ ],T in case C < 0. More precisely, using the standard 
symbol o(l) to signify a function such that lim^oo o(l) = 0, we expect that the 
equation has fundamental solutions of the form e ±Iv/ ^ T (l + o(l)) if C > 0, or 

e ±V / l^1 T (i + o(l)) if C < 0. Actually, this is the case if and only if the function 
W(t) satisfies the following finiteness condition (e.g. ,[7|); 

/oo />oo 
\W-C\dT= \W - C\adt < oo. (82) 

In the present case (with m ^ 0), it is confirmed that this condition is satisfied 
with the choice C = m 2 /4Z 2 > and 

a(t)=l-^. (83) 

With this, 

T = J <r(t)dt = t- 2^1og£, (84) 

and 

a{t) = \ + 0{^). (85) 
Recalling Eq. Q , we thus have the following: 

Proposition 2 The generic, m ^ 0, mode equation for massless scalar field 
on the LRS vacuum Bianchi VIII solution, satisfying Eq. \74\) , has the following 
fundamental solutions: 

X£\t) = t-i e ±*ISI(*-2piogt)( 1 + ^ (86) 

We must consider the m = case, separately. In this case, it is possible to 
choose 

a{t) = - t (87) 

so that 

W{t)=k-\ + 0{- t ). (88) 

Therefore, as long as A ^ 1/4, we can again apply the criterion mentioned 
above, since 

|W-(A-i)|a = 0(_L), (89 ) 

implying the integral 1)82(1 is finite with C = A — 1/4^0. Note that there are 
both possibilities of C being positive and negative, depending on which we have 
two cases: 

Proposition 3 The U (1) -symmetric, m — 0, mode equation for massless scalar 
field on the LRS vacuum Bianchi VIII solution, satisfying Eq. {7$ , has the 
following fundamental solutions: 

y(±wMer"M*«(l + (l)) (A>±) 

X ° {t) -\t-^(l + o(l)) (0<A<i), (90) 



where \s\ = y/\A - 1/4|. 
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We are interested in comparing with the Bianchi III model G2 • To 
compare the asymptotic solutions, it is convenient to use a proper time r as 
canonical time. Since our time coordinate t asymptotically approaches proper 
time, it might be expected that the formulas remain the same, but this is not 
exactly the case. As confirmed in the following, the phase velocity of a part 
suffers a slight modification. 

Let t be the proper time for the LRS solution, defined by Ea. i|tj(J|) . Using 
the formula (|61|l . it is easy to confirm, e.g., 

r i =T -i + (^), (91) 



and 



2/i log t = t + /x log t + 0(— ) - 2 M logT(l + 0(— )) 

T T 

= r-/ilogr + 0(^). 



(92) 



Beware that in the last equation the coefficient of log r is half the one of log t 
in the left hand side. Using these formulas, it is easy to rewrite the asymptotic 
solutions in Propositions [2] and [3] in terms of proper time r. As a result, we 
have the following. 

Theorem 4 A massless scalar field on a spatially compactified LRS vacuum 
Bianchi VIII solution can be decomposed into its mode components that are 
independent from each other and each of which is specified by the fiber eigenvalue 
m and the Casimir eigenvalue A. Their fundamental solutions (at the future 
asymptotics) are given in terms of proper time r as follows: 

r-5e ±i l s l 1 °s T (l + o(l)) (m = 0,A>l/4) (93) 

r -!±l s l(l + o(l)) (m = 0,0 < A < 1/4) 

Here, I and /j, are the parameters in the LRS solution 15H\I . and \ s\ = \/\A — 1/4 1 . 



6 Comparison with the Bianchi III model 

We make a comparison between the Bianchi VIII and Bianchi III systems. We 
first compare the backgrounds, and then compare the asymptotic behaviors of 
scalar field. 

The Bianchi III background solution. The LRS vacuum metric for 
Bianchi III is given by 

gpn) = -N\t)dt 2 + q^ma 1 ) 2 + (a 2 ) 2 ) + q 3 (t)(a 3 ) 2 (94) 

where 

N 2 {t) = t —^, qi (t) = (t-rf, g 3 (*)=4/ 2 ^, (95) 

and ii and I are real parameters of solution. The invariant 1-forms a 1 (I = 1 ~ 3) 
here are those for Bianchi III and must satisfy the following canonical relations; 

da 1 = a 1 A a 2 , da 2 = 0, da 3 = 0. (96) 
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Our convention in terms of coordinates is 



o~ = dx/y, a 2 = dy/y, cr 3 = dz. 
The invariant (dual) vectors are given by 

Xi=yd/ dx, X2 = y d/dy, xz = d/dz 



(97) 



(98) 



In this section vectors £/, xi> one-forms a 1 , and metric functions N(t), qi(t) all 
refer to those for Bianchi III (not for Bianchi VIII), unless otherwise stated. 

See Appendix A, , for details of the compactification of this solution. We 
assume an orthogonal |26| compactification, i.e., each hyperbolic z — constant 
surface is compactified to a higher genus surface S s . Although we do not repeat 
the details of the compactification, one of the most important facts for our 
discussion is that when compactified, each z-axis descends to circle fibers, and 
therefore the invariant vector X3 is called the fiber generator, like in the Bianchi 
VIII case. 

As easily confirmed from the above metric, the conformal base metric /r 111 ' = 
(cr 1 ) 2 + (cr 2 ) 2 is the same as the one for the Bianchi VIII; they both coincide 
with the standard hyperbolic metric h = (dx 2 + dy 2 )/y 2 . Since the 1-form cr 3 is 
dual to the fiber generator X3, Q3{t) is the scale factor for the fibers, while qi(t) 
is the one for the base. This interpretation of the metric functions qi(t) is the 
same as that for the LRS Bianchi VIII solution (j58(l . 

The above metric is essentially the same one adopted in 26 . However, we 
have made two alterations. One is that we have replaced the parameter k with 
—fi. As we will see, this makes the correspondence to the Bianchi VIII metric 
(|58|l better. The other alteration is the introduction of a new parameter I. At 
first, this parameter might seem redundant, since we could set 4^ 2 (cr 3 ) 2 — > (cr 3 ) 2 
by the induced map of the scaling diffeomorphism z — > z/(2l). However, we 
should notice that this diffeomorphism is not well defined on the compactified 
manifold. To see this, let us focus on the fiber submanifold (R, (^(cr 3 ) 2 ). To 
compactify we use the action by the translation z — » z + 2n. This map generates 
a group, 27rZ. The circle fiber T can therefore be expressed as 2-7rZ\(R, q 3 (a 3 ) 2 ). 
The translation however does not commute with the scaling diffeomorphism, 
which means that the scaling is not well defined on the compactified manifold. 
The parameter I is therefore not redundant for the compactified manifold. The 
significance of this parameter is apparent if we consider the length of the fiber, 
which is given by 



In particular, Loo = limbec L^ 111 ' = Anl. Therefore I is the (47r) _1 times the 
fiber length at the future infinity, like the LRS Bianchi VIII case. (See Ea . (1631) . ~) 
We comment that we could also use the fiber metric q' 3 (a 3 ) 2 that is obtained 
by setting 4/ 2 = 1 as the universal cover metric, to express the compactified 
fiber. To this, we need to allow the covering group (rather than the metric) to 
have the parameter I, i.e., to compactify the fiber we consider the one-parameter 
translation z — > z + 4ttZ. Then, it is easy to see that the resulting fiber T' = 
4irlZ\(R, <?3(cr 3 ) 2 ) is equivalent to the fiber T . For example, the fiber length is 
computed as 




(99) 




(100) 
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which agrees with L( ln \ Remark that the parameter I is a relevant parameter in 
any case. In general, given a compactified locally homogeneous manifold, there 
are two ways to express it depending upon whether we fix the covering group or 
not. See [2HJ for a treatment of spatially compactified models with fixed covering 
groups (i.e., with no dynamical degrees of freedom in the covering group acting 
on spacetime), and )24j for one with varying covering groups with the smallest 
number of parameters in the universal cover metric. (See also 1131 112) for 
related discussions.) Note that in the background spacetime is expressed 
in the latter view point, while in this paper the former point of view has been 
exploited. 

The reason we introduce the two-parameter metric (19511 is that it has a 
straightforward correspondence to the (LRS) Bianchi VIII metric (|58f) . In par- 
ticular, as mentioned, the parameter I has the same meaning that it is (47r) _1 
times the fiber length at the future infinity. 

Let r be the proper time defined by 



w< Wv^W( i -i +0( ^ ,|<a 

t- fi\ogt + 0{- t ). 



(101) 



This implies 



t = r + /ilogr + 0(^). (102) 



Therefore, from Eos. 1)95)1 we have 



= (r + Mlogr) 2 + 0(logr), 

, 3(T) = W 1 + ^ + °(^>. (103) 



T J 



T 2 



Comparing these equations with Eas. 1)62)1 it is confirmed that the asymptotic 
behavior of the LRS Bianchi III model is the same as that of the Bianchi VIII 
at least up to second leading terms. That is, the future asymptotic behaviors 
of base surface and circle fiber are the same for both LRS Bianchi III and VIII 
solutions. This establishes the background degeneracy of the two models. 

Asymptotics of scalar fields on the Bianchi III. To compare the be- 
haviors of scalar field we need to summarize the asymptotic properties of scalar 
field on the LRS Bianchi III background. As far as the compactified spatial 
manifold is of the orthogonal type, the harmonics can be constructed by simply 
making products of those c m (z) on the fiber and those S\(x,y) on the base 
hyperbolic surface. Those harmonics are defined by the following eigenvalue 
equations 



X3C m = imc m , 

2" ( 104 ) 

\S\ = —A S\. 



Here, A is the Laplacian with respect to the standard hyperbolic metric h — 
(cr 1 ) 2 + (a 2 ) 2 . The harmonics (mode functions) on an orthogonal closed Bianchi 
III manifold are the products 

S\, m = c m S\. (105) 
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We call m the fiber eigenvalue, and A 2 the base eigenvalue. Remember that 
since we take the two-parameter metric (|95|) . the identifications along the fibers 
are taken with the fixed step 2n, and as a result we have 

meZ, (106) 

since the solution of the first equation in Eas. p()4(l is given by e lmz . 

The mode-decomposed massless scalar field equation, computed with the 
one-parameter metric, is given in Eq.(243), |26|. To convert this equation to the 
one computed with the two-parameter metric l|95|) . it is enough to perform the 
replacement 

Then, the equation reads 

^^+G^ + (i) J trf> = °- « 

The unknown function ip = ip\ tm (t) is the field component for the mode 5a, m- 
The exact solutions for the /i = background are given in Appendix B, |26| . 
No results are however presented there for general [i. Fortunately, however, it is 
not difficult to obtain asymptotic solutions for the equation, following the same 
procedure shown in the previous section. We just present the result here. 

Proposition 5 Consider the massless scalar field eguation on the LRS Bianchi 
III vacuum metric {9$ , given by Ea. l\lU8\ ). The fundamental solutions of the 
equation are given by 

Yj^Ht) = \ t-5e ± *l s l 1 °s*(l + o(l)) (m = 0,A 2 >l/4) (109) 

1^-5^1(1+0(1)) (m = 0,0< A 2 < 1/4), 

where 



|A 2 -j|. (HO) 

It is easy to convert t into the proper time r in this solution. Using Ea. H102|l . 
we obtain the following. 

Theorem 6 A massless scalar field on a spatially compactified LRS vacuum 
Bianchi III solution can be decomposed into its mode components that are in- 
dependent from each other and each of which is specified by the fiber eigenvalue 
m and the base Laplacian eigenvalue A 2 . Their fundamental solutions ( at the 
future asymptotics ) are given in terms of proper time r as follows: 

r -l e ±i|f ICr-Mlcgr)^ + ( TO ^ ) 

>1 ±) W = { r- = e ±i l s l 1 °s T (l + o(l)) (to = 0,A 2 >1/4) (HI) 

3±l s l(l + (l)) (to = 0, < A 2 < 1/4), 



Here, k and I are the parameters in the LRS solution 195\) . and \ s\ = ^/\\ 2 — 1/4|. 
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Again, only the difference from the version expressed in terms of the coordinate 
time t is that the numerical factor of — 2filogt, in the phase part of the solution 
for m 7^ 0, becomes half, — /ilogr. 

Comparison. Comparing the above theorem with Theorem we see that 
the asymptotic solutions in the Bianchi VIII and Bianchi III models completely 
agree, provided the correspondence A 2 <-> A is understood. We call this agree- 
ment the scalar field asymptotic degeneracy of the two models. 

In the following we give a short account how this degeneracy can be un- 
derstood. First, note that an arbitrary scalar field solution \& can be de- 
composed into two parts, the U(l)-symmetric part Y(t,x,y) and the rest part 
^( g en){t,x,y,z) = * - T: 

* = *(gcn)+T. (112) 

(We do not have to mode-decompose each part here.) We first consider the U(l)- 
symmetric part T. Remember that the U(l)-symmetry means the translation 
symmetry along the fibers. As mentioned, as long as the background is LRS 
the background is also U(l)-symmetric along the fibers for both Bianchi types. 
Therefore we can consistently contract the fibers and obtain a reduced scalar 
field system on a (2 4- l)-dimensional spacetime. The 2-dimensional spatial 
manifold is the higher genus surface E a , and the spacetime manifold becomes 
R x E s for both models. It is easy to write down the field equation on this 
contracted manifold for the Bianchi VIII system. From Eas. (|65|) . i|titj|) and (|53[) . 
we have 

g afc V a V b T = (-^L| (V^V- 2 |) + ft 1 A ) T, (113) 

where A is the Laplacian (|55|) with respect to the standard hyperbolic metric. 
Thus, the field equation on the contracted manifold is given by 

N d ( qiy/qj d 



— . , , , ■ N V An T = 0. (114) 

In this equation, the functions qi(t) and N(t) should be thought of as the func- 
tions of time given in Eq. (|58|l . The equation for the Bianchi III is also obtained 
in the same way, and we find that the above equation is exactly valid for this 
case, too, provided that qi(t) and N(t) are those for the LRS Bianchi III so- 
lution i|95|) . This equivalence explains the degeneracy for the U(l)-symmetric 
modes, since we have the background degeneracy, which means that the coef- 
ficient functions in Ea. H114(l show asymptotically the same behavior for both 
cases. 

Next, consider the rest part ^(gen), which we call the "generic part." Note 
that the spatial Laplacians for the two LRS backgrounds are expressed as 

a _/?r 1 ((xi) 2 + (x2) 2 ) + ^ 1 (x3) 2 (viii) 
9 " Vao + ^W. (in) 



The key fact is that in bath cases, q 3 dominates q 1 , J 

l im «iW = o. (116) 
% \t) 

5 It is noteworthy that the dominance of q^ 1 also holds for the non-LRS Bianchi VIII 
background, i.e., as we can check using a result in |17) . it holds limq^ 1 /q^ 1 = limij^" 1 /q^ 1 = 
0. 
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Since from the assumption we have X3^(gen) 7^ 0> this may suggest that the 
term in q^ 1 will dominate the term in q^ 1 in both Laplacians. If this can be 
justified, we have the simplification 

A q ^q^( X 3) 2 , (VIII and III) (117) 
resulting in the same asymptotic equation for both systems again; 

(-^(^) +w2 ^*->^ <118) 

This explains the degeneracy for the generic modes. 

We call the last equation the fiber term dominated (FTD) equation of the 
scalar field equation. 6 Beware that in contrast to the U(l)-symmetric equation 
(|114Jl . the above equation cannot be justified as a symmetry reduction. That 
is, although we dropped off \i an d X2~dependent terms, it is not to impose the 
additional conditions xi^ = X2^ = 0. Indeed, if this is the case for the Bianchi 
VIII we must have = — [XijX^]^ = 0, implying only spatially constant 
configurations are allowed. Therefore \& has no symmetry and in particular 
depends on all spatial coordinates "J" = ^(t, x, y, z). 

Note that although Ea. (|117|l looks natural, we need a proof for its justifica- 
tion, which is what we have done with the mode decomposed equations. 

7 Conclusions 

We have separated the scalar field equation on a Bianchi VIII background using 
the harmonics for the circle-fibered closed Bianchi VIII manifold. The reduced 
wave equations form two sets of infinite number of simultaneous ordinary dif- 
ferential equations (ODEs) . This is a result of the fact that irreducible repre- 
sentations of a noncompact Lie group like the Bianchi VIII group are in gen- 
eral infinite dimensional. When the background is LRS, however, each single 
wave equation becomes closed itself due to the additional symmetry. We have 
analyzed this closed equation and obtained future asymptotic solutions, as in 
Theorem^] In particular, for the m = (U(l)-symmetric) and A < 1/4 ("small 
eigenvalue" ) case the solution is non-oscillatory, while the other cases are oscil- 
latory. They are all decaying (except one of the two fundamental solutions for 
the "zero- mode" with m = A = 0). 

We have seen that this result completely agree with that of the Bianchi III 
model. We interpret this "scalar field asymptotic degeneracy" as an evidence 
that other linear fields, including electromagnetic fields and linear perturba- 
tions, on the Bianchi VIII background also have the same asymptotic behaviors 
as those on the Bianchi III. Since many results are already known for the Bianchi 
III we can i based on them, conjecture corresponding properties. Of our 

interest is that of the linear perturbations. The asymptotic solution of the sta- 
bility measures (a kind of normalized gauge-invariant variables) for the Bianchi 
III, given in Theorem 2.5, [22], shows that the Bianchi III vacuum solution is 
asymptotically unstable, meaning that the perturbed spacetime asymptotically 

6 This equation can be considered as an analogy of the asymptotically velocity term domi- 
nated (AVTD) equations 1101 of Einstein's equation valid for the opposite time direction. 
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becomes more and more inhomogeneous. It is therefore plausible that the (LRS) 
Bianchi VIII vacuum solution is also asymptotically unstable in the same sense. 

One however needs to be careful about the difference between the properties 
of the vector or tensor harmonics applied to the two models. For example, 
the Bianchi III tensor harmonics arc split into four kinds, the "even" ones, the 
"odd" ones, the "harmonic" ones, and the "transverse-traceless (TT)" ones |2fci| . 
Accordingly, there are four kinds of (independent) perturbations in the Bianchi 
III case. On the other hand, we will have no such splitting in the Bianchi VIII 
case. This difference comes from whether the invariant frame {a 1 ,xi} is wen 
defined on the compactified manifold considered, since if it is well defined we 
can construct all vector or tensor harmonics from the scalar harmonics with the 
help of the frame. See [23 for an explicit example, where a Bianchi II case is 
treated. Because of this property, the Bianchi VIII tensor harmonics will not 
have a splitting like the Bianchi III one. 

Because of this difference, the reduced Bianchi VIII perturbation equations 
will have somewhat different properties. In particular, while, in the Bianchi III 
case, the reduced perturbation equations are given as independent second order 
ODEs (or independent systems of two simultaneous first order ODEs), in the 
Bianchi VIII case they will be given as independent systems of two simultaneous 
second order ODEs (or independent systems of four simultaneous first order 
ODEs). (Here, we are assuming an LRS background for both Bianchi types.) 
We can interpret this increase of variables in an independent system of equations 
as a consequence of the coupling between an "even" and "odd" mode, like the 
Bianchi II case |23| . 

The result we should refer to to obtain a corresponding property about the 
Bianchi VIII perturbations should be the one for the "even" perturbations of 
the Bianchi III solution. This is because they are the dominant perturbations 
among the even and odd ones. It is apparent that the Bianchi VIII perturbations 
will have nothing to do with the "harmonic" and "TT" ones of the Bianchi III 
solution, since those perturbations are connected to the modes that cannot be 
produced from the scalar harmonics. As a result, we can conclude that the 
conjectured growth rate of the stability measure for the generic Bianchi VIII 
perturbations is 0(t) (Cf. Theorem 2.5, [22]), meaning it is unstable. 

Finally, we comment on non-LRS cases. Remember that the degeneracy 
between the Bianchi VIII and III systems concerns the LRS backgrounds. Al- 
though, since a Bianchi III manifold cannot be compactified unless it is LRS, we 
do not consider non-LRS Bianchi III system, the non-LRS Bianchi VIII system 
may be of great interest itself. However, the apparent difficulty is that we must 
solve infinite number of simultaneous ODEs in this case. 

We point out that this difficulty might not be so crucial. We have seen 
that the scalar field asymptotic degeneracy (for the generic modes) was a con- 
sequence of the fiber term dominated (FTD) behavior. Since the dominance of 
the scale factor function q^ 1 over q^ 1 and q^ 1 continues to hold for the non-LRS 
background, we can expect that an FTD behavior holds for this case, too. If this 
is justifiable, it implies that each single mode becomes (virtually) independent 
like in the LRS case. (This is what can be easily confirmed by inspecting the 
non-LRS case equation.) Therefore the FTD behavior may be the key to find 
asymptotic solutions of the non-LRS case field equations. 
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